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One of the most challenging issues in the characterization of magnetic materials is to obtain 
quantitative analysis on the nanometer scale. Here we describe how electron magnetic circular 
dichroism (EMCD) measurements using the transmission electron microscope (TEM) can be used 
for that purpose, utilizing reciprocal space maps. Applying the EMCD sum rules, an orbital to spin 
moment ratio of rriL/ms = 0.08 ± 0.01 is obtained for Fe, which is consistent with the commonly 
accepted value. Hence, we establish EMCD as a quantitative element specific technique for magnetic 
studies, using a widely available instrument with superior spatial resolution. 
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Fast advances in the field of magnetic nanostructures, 
both in fundamental research and technological develop- 
ment, call for new magnetic characterization methods. 
Electron microscopy is nowadays a standard technique 
for structural and chemical analysis down to the atomic 
scale. Magnetic imaging 1 in the TEM is also possi- 
ble while measurements of element-specific magnetic mo- 
ments have until now been the domain of synchrotron 
based dichroic experiments, such as x-ray magnetic cir- 
cular dichroism (XMCD) [2]. Although XMCD is widely 
applied in materials science, it is mainly related to surface 
measurements and with limitations in spatial resolution. 
The work of Schattschneider et al. [2j - reporting an ob- 
servation of dichroic effects in the TEM ~ opened a new 
route for high-resolution element specific magnetic char- 
acterization, using widely accessible standard laboratory 
equipment . 

EMCD measurements are in principle simple. An un- 
polarized electron beam, passing through a magnetic ma- 
terial, exhibits a magnetic dichroism in the momentum 
resolved electron energy- loss spectra (EELS) |3]- The 
origin of this effect stems from the inelastic scattering of 
incoming high-energy electrons that excite core electrons 
to unoccupied states. The signal at a scattering vector 
k contains mixed contributions of all pairs of diffracted 
beams with momentum transfers q and q'. A dichroic ef- 
fect appears when two EELS-spectra - extracted at spe- 
cific detector positions in reciprocal space defined by a 
mirror axis - are subtracted (see Fig.[l|. This difference 
spectrum is called in the following the EMCD signal. 

While the principle of EMCD has been demonstrated, 
decisive progress is required to allow quantitative mag- 
netic analysis, which is reported here. The recent deriva- 
tion of the EMCD sum rules for extraction of spin (mg) 
and orbital (wi) magnetic moments represents an im- 
portant step in that direction [U [3- As EMCD relies 



on reciprocal space vectors, proper k-space selection of 
detector positions is essential. So far, most measure- 
ments are carried out by selecting a limited part of recip- 
rocal space from where the EELS-spectra are acquired 
[31 [6]. As shown in this letter, increased flexibility for 
data optimization and precision in k-space selection is 
obtained when a map showing the distribution of the 
signal in reciprocal space is available. Energy filtered 
reciprocal space maps can so far only be considered as 
semi-quantitative |5j|7j, since only raw-data was consid- 
ered. For a quantitative measurement, the signal to noise 
(S/N) ratio must be substantially increased and a set of 
procedures for data analysis of core-loss edge energy in- 
tensities must be devised. 

We choose to demonstrate the technique on bcc Fe, pri- 
marily because its magnetic properties are well known, 
which allows for a precise exploration of the EMCD tech- 
nique. Therefore, a Fe sample grown by UHV mag- 
netron sputtering (texture angle of ±0.3°) was chosen 
for the analysis. Three-dimensional data sets, so called 
data cubes, consisting of the reciprocal k^ — ky plane 
and electron energy-loss are acquired to obtain maps of 
the EMCD signal. The data cubes are similar to the 
spectrum-images introduced by Jeanguillaume and Col- 
liex [5], but here diffraction patterns instead of images 
are used. The energy filtered diffraction patterns were 
acquired using a Gatan GIF2002 spectrometer on a FEI 
Tecnai F-30ST FEG operated at SOOkV with energy step 
of leV (slit width 2eV) at a sample region of thickness 
19 ± 2nm and an electron beam diameter of 600nm. 

The geometrical conditions of the EMCD experiment 
are crucial for understanding the effect. The first EMCD 
experiments [5] were performed with the sample oriented 
in a two beam case (2BC) geometry, see Fig.fll The sam- 
ple is tilted with respect to the incoming electron beam 
from a high symmetry orientation by an angle a ^ 10°, 




FIG. 1: A sketch of the diffraction geometry and distribution 
of EMCD signal in the k^ — ky plane. The sample is tilted 
with respect to the incoming electron beam from a high sym- 
metry orientation to the 2BC geometry where a ~ 10° and 
/3 « 0.4°. In this geometry only the transmitted (0) and 
Bragg scattered beams (G) are strongly excited (large black 
spots). The optimum detector positions are indicated as PupR 
and PonR for 2BC geometry using a horizontal mirror axis 
(solid line). For the SBC geometry, a « 10° and /3 = 0°, an 
additional vertical mirror axis is available (dashed line). In 
the experiment, energy filtered diffraction patterns appear in 
the kx — ky plane; here an EMCD map is shown as a guide to 
the eye. The black circle indicates detector positions where 
q _L q', i.e., the Thales circle. 



where mainly a row of reflections in reciprocal space is ex- 
cited. By tilting the sample further in the perpendicular 
direction by a small angle oi (3 ^ 0.4° the 2BC geom- 
etry is obtained. In this geometry the transmitted and 
Bragg scattered beam G = (200) in Fe are strongly ex- 
cited, while all others are weak. To use the sum rules [1], 
allowing a quantitative assessment of the mLJms ratio, 
the dichroic signal must be extracted using two symmet- 
rically placed detectors in reciprocal space and the in- 
coming beam should lie within a mirror symmetry plane 
of the crystal structure. In 2BC geometry the EMCD 
signal is given by the difference of a spectrum extracted 
in the upper half plane (at PupR) and the corresponding 
spectrum in the lower half plane (at Pdur) in recipro- 
cal space, see Fig. [T] This defines a mirror axis for the 
2BC geometry, denoted as horizontal mirror axis. By 
using this mirror axis and correspondingly subtracting 
all spectra in reciprocal space, maps of the EMCD sig- 
nal are constructed. However, the 2BC geometry does 
not fulfill the condition of symmetric detector positions. 
This can be seen in Fig. IT] where the (200) atomic planes 
are not symmetric with respect to the positions of the 
detectors {PjjpR and Pdur)- Then the cancellation of 
non-magnetic contributions to the signal is not perfect 



[HI HO] , leading to limitations in the use of the 2BC ge- 
ometry for mhlfns ratio determination. Therefore, we 
apply the three beam case (SBC) geometry, which is fully 
symmetric. Here, the incoming beam is oriented within 
a symmetry plane (200), /3 = 0° in Fig.fll The condition 
of symmetric detector positions is now fulfilled when ex- 
tracting the dichroic signal as a difference between left 
and right half plane, defining a vertical mirror axis. The 
SBC geometry also enables the use of a horizontal mirror 
axis. Using both mirror axes improves the S/N ratio by 
using data from all four quadrants. Moreover, it can be 
shown that such a double difference procedure corrects 
slight misorientations from exact SBC geometry |10j . 

The EMCD signal evolves as an entwined property of 
sample thickness and orientation due to dynamic elec- 
tron diffraction effects. Thus, simulations of the EMCD 
experiments include both Bloch-wave simulations of dy- 
namic scattering of the fast electrons, as well as ab ini- 
tio calculations of the inelastic scattering in a magnetic 
sample |llj . Inelastic transition matrix elements were 
calculated including spin-orbit coupling with magnetiza- 
tion direction [016] as imposed by the magnetic field of 
the objective lens, corresponding to a « 10°. In order 
to simulate the finite collection angle and the relative 
misorientation of illuminated regions in the experiments 
(±0.S°), maps of the EMCD signal were averaged over 
a set of Laue circle centers. Simulated maps of the rela- 
tive EMCD signal at L3 edge for 2BC and SBC geome- 
tries are displayed in Fig. [2^ and [2^, respectively. Since 
DFT calculations predict a very low orbital momentum 
(~ 0.045^b) in bcc Fe, the EMCD map at L2 has a 
similar structure as the L3 map but with reversed sign 
PT| . In the 2BC geometry, the strongest dichroic signal 
is shifted towards the strongly excited refiection, denoted 
as G in Fig. Ill in agreement with Refs. [12l [13]. A re- 
duced EMCD signal is also present around the weakly 
excited reflection (— G). 

After applying cross-correlation [T?] on the transmit- 
ted beam, each spectrum in the experimental data-cube 
was pre-edge background subtracted (power-law model) 
and normalized in the post-edge region (at 7S6-7S8eV 
as seen in Fig. Isk). After that, each spectrum in the 
kx ~ ky plane was peak fitted |15j and the integrated L2_3 
edge intensities are used to construct the EMCD maps. 
Peak fitting removes the need of somewhat arbitrary en- 
ergy window selection for separating L3 and L2 compo- 
nents and corrects for eventual energy shifts due to non- 
isochromaticity of the spectrometer. Experimental maps 
of the relative EMCD signal at the ^2,3 edges in 2BC 
and SBC geometries are shown in panels b, c, f and g of 
Fig. |2] respectively. We find a very good correspondence 
between the theoretical prediction (panels a and e) and 
measurements regarding the position in reciprocal space 
and strength of the EMCD signal for both geometries. 
From the 2BC data cube two spectra were extracted us- 
ing a box with the size of O.5G200 x O.5G200 at the PjjpR 



and Ponfl-positions (similarly to Ref. [T3]) and peak fit- 
ted. In Fig. Isk the ^2,3 edge intensities and a continuum 
background used in the fit of the experimental data are 
shown. The spectra show a clear EMCD signal with good 
S/N ratio. Experiment and simulation reveal a relative 
EMCD signal (difference divided by its sum) at the L3 
edge of approx. 8% and 12%, respectively. Note that a 
strong EMCD signal is obtained in spite of the texture 
angle of ±0.3° = ±5mrad (corresponding to ±0.4G2oo), 
which is of high importance for a practical use of EMCD. 
Thus the technique is robust against small distortions of 
the crystal orientation, allowing for the characterization 
of imperfect crystals. A texture angle will have a simi- 
lar effect on the EMCD signal as an increase of conver- 
gence angle since both can be viewed as electrons hitting 
a single crystalline lattice at slightly different angles. A 
semi-convergence angle of 10 mrad corresponds in the 
TEM used here to a probe size of about 0.13 nm [16]. In 
this work, wc show that a quantitative EMCD signal is 
obtainable at deviations in the relative sample to probe 
orientation of ±5mrad. Therefore we believe that the 
approach demonstrated here has the potential for sub- 
nanometer resolution. This argumentation is underlined 
by the recent demonstration of an EMCD signal using 
a convergent electron beam at 2nm resolution |17j . The 
quest to reach the ultimate resolution of this technique 
will thus not be dominated by the achievable convergence 
angle or electron probe size, but by the stability of the 
sample and the electron microscope itself. 

A quantitative assessment of spin and orbital angu- 
lar momenta can be obtained from the EMCD sum rules 
[31 |S]. Correspondingly, a map of the mi^jms ratio is ob- 
tained by calculating at each position in reciprocal space: 

ms 3 L M{E)dE - 2 /^, M{E)dE 



(1) 



with AI{E) = I{Pi) - I{Pj) where I{Pi) and I{P,) are 
the EELS-spectra extracted at the two detector positions, 
e.g. in the 2BC geometry a choice of Pi = PupR imphes 
using Pj = PduR see Fig. IT] The general sum rule ex- 
pressions [4 predict a k-independent mL/ms ratio, i.e. 
the intensity of a large part of reciprocal space can be 
employed to determine the ratio. Nevertheless, the S/N 
ratio should be optimized by selecting the region with 
highest EMCD signal. Since sum rules involve several 
approximations and assumptions, which may not be per- 
fectly fulfilled, we evaluated theoretical tol/tos maps 
(not shown). They reveal the effect of 2BC asymme- 
try [9 , which becomes non-negligible in regions outside 
the Thales circle. In the 3BC geometry even a devia- 
tion as small as /3 = O.O5G200 can introduce substantial 
variations of the observed mL/ms ratio throughout the 
diffraction plane. Detailed calculations have shown that 
construction of the double difference map, i.e., exploiting 
both mirror axes, is a very efficient method for correcting 
these inaccuracies [TO] . 




FIG. 2: Reciprocal space maps of the EMCD signal and 
rriL/ms ratio for the Fe sample oriented in a 2BC geometry 
(left) and SBC geometry (right). Theoretical relative EMCD 
maps at the L3 edge are shown in a) and e). The inset shows 
the simulated diffraction pattern. In b) and c) [f) and g) for 
3BC] maps of experimentally obtained relative EMCD sig- 
nal at L3 and L2 edges are shown. The black lines indicate 
the applied mirror axes and blue spots the positions of the 
transmitted and Bragg scattered G = (200) and — G = (200) 
beams. The insets in b) and f) show the diffraction patterns 
averaged over an energy interval from 695 eV to 740 eV. In d) 
and h) the experimental mL/ms maps are shown, where the 
box indicates the region with minimal noise. Values outside 
the range of color bar are shown as white. 



The maps enable an accurate determination of the 
'm,L/'ms ratio. A selection window is used in these maps 
from which the values of all individual pixels are ex- 
tracted, thus optimizing the S/N ratio. The overall 
'niL/'ms ratio for a given window size (collection angles) 
is obtained by fitting the histogram of the individual 
TnL/ms values with a Gaussian (see inset in Fig. pp). 
The ratio is determined as the center of the Gaussian fit. 
Plotting this ratio as a function of the window size, a close 
to constant value in the m-L/rng ratio is reached when the 
size of the box is sufficient to reduce the statistical fluc- 
tuations. If the window gets too large, regions with low 
S/N ratio are included which increase the standard de- 
viation of the rfiL/rnig ratio. Regions of minimal noise 
in the experimental vfiL/rng maps (indicated in panels d 
and h of Fig. p| correlate with local maxima found in the 
EMCD maps. These results agree well with findings in 
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FIG. 3: Experimental EELS-spectra at Fe L2,3 edges and 
rriL/ms ratios as a function of collection window size in 
'/nL/ms ratio maps. The two spectra in a) were extracted 
from the experimental data cube with the sample in 2BC 
geometry at the PonR (orange squares) and PupR positions 
(blue circles), averaging over O.5G200 x O.5G200 in reciprocal 
space. The components of the fit function for the spectrum at 
the PupR position are shown, b) mL/ms ratio obtained for 
various window sizes in 2BC geometry using horizontal mirror 
axis and in SBC using double difference method. Inset shows 
the histogram with fit of rriL/ms ratio for the window size of 
O.5G200 X O.5G200 in 2BC orientation (indicated in Fig.[2}i). 



Ref. [ISj , although here we apply the integration window 
numerically on the mi^/mg ratio rather than on maps 
of the EMCD signal. This provides better flexibility for 
numerical analysis and a reliable error bar. Here, an op- 
timum size of the window in reciprocal space is typically 
0.56*200 X 0.56*200- We obtain a consistent rriL/ms ra- 
tio, depicted in Fig. |3]d, of 0.09 ± 0.01 in the 2BC and 
0.08 ± 0.01 in SBC geometry using the double difference 
maps. The standard error s — 0.01 was estimated using 
N = 1225 individual rriL/ms ratios within the selection 
window (standard deviation of an individual rriL/ms ra- 
tio is cr=0.23, s=a/VN). 

Our experimental rriL/rns ratio for bcc Fe is 0.08 ± 
0.01. The most commonly accepted values of the ratio 
are obtained by XMCD and gyromagnetic ratio measure- 
ments, 0.043 IISJ and 0.044 [19], respectively. One may 
observe that the present value is somewhat larger. How- 
ever, neutron scattering experiments give a value closer 
to the presently obtained value, namely 0.062 [2D]. For 



bcc Fe the orbital moment is small and it is not surpris- 
ing that this delicate property shows certain dispersion. 
E.g., with XMCD, values of 0.043 [18], 0.07 [21], 0.085 
[H] and 0.12 ± 0.05 [55] for the mi^/ms ratio have been 
reported for bcc Fe. In the experiments presented here, 
oxygen is most likely present in the vicinity of one of the 
Fe surfaces, which may also influence the tol/tos ratio. 
It should be noted that one previous EMCD experiment 
measuring the mi^/mg ratio has been conducted for Fe 
[S]. Although the rriL/ms ratio reported in that work 
is close to ours, their analysis lacked a quantitative data 
treatment (e.g., peak fitting), as they stated explicitly. 

Experimental and theoretical maps of the EMCD sig- 
nal in reciprocal space are shown to agree well for the 
investigated sample geometries. Asymmetry problems 
are successfully compensated by the double difference 
method within the 3BC geometry. Obtained mL/ms ra- 
tios of bcc Fe are close to the commonly accepted value. 
Our work establishes therefore the quantitative use of 
EMCD, which, together with the high spatial resolution 
available in modern TEMs promises unique possibilities 
for truly nano-scale magnetic studies. 
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